Nonlinear Wave Interactions as Emission Process of Type II Radio Bursts 
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ABSTRACT 

The emission of fundamental and harmonic frequency radio waves of type II radio 
bursts are assumed to be products of three-wave interaction processes of beam-excited 
Langmuir waves. Using a particle-in-cell code, we have performed simulations of the 
assumed emission region, a CME foreshock with two counterstreaming electron beams. 
Analysis of wavemodes within the simulation shows self-consistent excitation of beam 
driven modes, which yield interaction products at both fundamental and harmonic 
emission frequencies. Through variation of the beam strength, we have investigated the 
dependence of energy transfer into electrostatic and electromagnetic modes, confirming 
the quadratic dependence of electromagnetic emission on electron beam strength. 

Subject headings: plasmas, Sun: heliosphere, Sun: radio radiation 



1. Overview 



Emission of radio b ursts from the solar corona has been observed since the early 1950ies 
(jWild &: McCreadvl Il950l ) . While theoretical explanations for the emission mechanisms of most 
types of bursts have been found, many details of the emission process of type II radio bursts are 
still a matter of discussion. 



The morphology of type II bursts typically shows a two-band emission spectrum, consisting of 
the fundamental emission band and of the harmonic emission band at about twice the frequency 
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of the fu ndamental (in a few , near- limb events signals of third harmonic emission are also discern- 
able, see lZlotnik et al.lll998l ). The fundamental frequency is believed to correspond to the plasma 
frequency of the emission region, slowl y decreasing over time as the coronal/int erplanetary shock 
travels outwards into the heliosphere. (jCane et al.lll987l ; iNelson Melrose! 1 19851 ) 



The commonly accepted model by iReiner et al.1 (119981 ) proposes that the emission region is 
located around points where the CME drives a curved perpendicular shock, leading to efficient 
shock drift acceleration of electrons. Due to the curved nature of the shock, these can escape 
into the foreshock following the magnetic field, whe re they form an el ectron beam population. 
This model has been thoroughly treated analytically. iKnock et al.l (|200ll ) have derived theoretical 
emi ssivities for Langmuir wave emi ssion and subsequent conversion into electromagnetic emissions, 
and lSchmidt Gopalswamyl (120081 ) have combined these with MHD simulations to track the radio 
burst emission region as it moves outward from the Sun. These theoretical predictions are however 
not in complete agreement with observations, and many details of the fine structure ( band splitting 



and herringbone patterns) often visible on type II bursts are yet to be explained (lAur ass et al 
1994T ). 



Similar to the assumed emission process of type III bursts, instabilities driven by the electron 
beams can then lead to creation of Langmuir waves (L) and their three wav e interaction products 
(Ion-soundwaves S and transverse electromagnetic waves T) (lMelrosdll986l ): 



L -± L' + S 
L S + T(u pe ) 

L + L' -»■ T(2co pe ) 



(1) 
(2) 
(3) 



Momentum an d energy conservation as well as wave polarization impose certain limits on the 
waves' k directions (jSpanier fc Vainidl2009l ). The Langmuir waves (L) should primarily be excited 
parallel to the beam dire ction, potentially following the beam electrons' pitch angle distribution 
(jKarlicky h Vandasl 120071 ) . whereas electromagnetic waves are expected to be emitted in quasi- 
perpendicular directions: Assuming that 



jfej- = fcjf + fcf 

0J L = LO T +UJ S , 



(4) 
(5) 
(6) 



Angular momentum conservation leads to polarization selection rules (jSpanier fc Vainioll2009i ). 
Inserting the dispersion relations for all three waves and using these selection rules to neglect k^ 
as well as k^_, we get 



3A;||f th = c kj_ + k^v s + 2k\\v s uj pe + k±k\\v s c, (7) 
(with speed of sound v s and electron thermal speed vth) which results in an expected A;-value for 
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the harmonic electromagnetic emission of 



< 3 jfc2(^L_ JJ< 
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(8) 



Being a second-order coupling of Langmuir waves, the intensity of the electromagnetic emission is 
expected to be quadratic in the Langmuir wav e energy, and should thus also be quadratic in the 
initial electron beam energy (jKnock et al.ll200ll ) for the model described here. This assumes a fixed 
level of ion soundwave intensity, which the simulations in this work provide through identical ther- 
mal initialization of the background ion population in all runs. Excitation of additional soundwave 
intensity by the beam can be neglected due to the small simulation timescales. 



In the simplest case of this emission scenario, a single electron beam emission site, created on a 
smoothly curved CME shock, is supplying the electron beam populations for radio emission. Lang- 
muir waves excited through the beam driven instability would first have to scatter off soundwaves 
(Eq. [[]) at a sufficient rate to create a suitable population of counter-propagating Langmuir waves 
L'. If this process is too slow, harmonic emission (Eq. [3j) cannot be expected to attain intensities 
which are comparable to the fundamental plasma emission process. 

An alternative scenario can be constructed by assuming not one single, but two (or multiple) 
points along a shock front in which perpendicular shock conditions pertain, allowing for drift 
acceleration of electrons into the foreshock. This could be caused by shock ripples at the front of the 
CME, or density fluctuations in the background solar wind leading to inhomogeneous propagation of 
the shock through the heliosphere. With multiple electron acceleration sites feeding electron beams 
into the foreshock from different directions, counterstreaming electron beams may be present within 
the emission region. Direct excitation of counter-propagating Langmuir waves through beam driven 
instability is hence possible, leading to a stronger harmonic emissivity. 



Due to the relatively small size of the emission region, in-situ satellite observation data is 
very scarce. In a handful of cases however, measurements of interplanetary type II burst emis- 
sion regions were possible , as an outward-expanding CME shock front encountered spacecraft at 
lAU ()Pulupa fc Baldl2008l ). The measurements obtained from these emission-region crossings agree 
with the proposed model, in that they provide evidence of counterstreaming foreshock electron 
beam presence. Electric field measurements on board the spacecraft are also indicating a strong 
presence of waves near the plasma frequency, attributed to Langmuir waves. However, direct in- 
situ satellite measurements of wave quantities are limited to the resolution of io s in the spacecraft's 
frame. S ince the solar wind is streaming by the spacecraft at a velocity Usw> 

this quantity is a 

mixture of to s = to + u sw ■ k, but the k contribution is negligible for the high-frequency wavemodes 
of interest to this study. While multispacecr aft obser vations have been successful in reconstructing 
directionalities of radio bursts (jTheiappa et al.ll2012l ). a detailed analysis of the kinematics of three 
wave interaction processes, for which information with both to and k resolution is essential, can 
only be provided by simulations. 
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1.1. Related Work 



In similar work. iKarlickv Vandaa (|2007l ) performed 1.5D simulations of beam-excitation of 
Langmuir waves from foreshock electrons, and subsequent couplings to kinetic wavemodes. Due to 
the limitation to one spatial dimension, these runs were not able to represent waves with transverse 
m omentum and hence prec luded the subsequent electromagnetic emission mechanism (eqs. ElEj). 



In lKarlicky &: Bartal (120101 ). 2.5D PiC simulations were conducted into which monochromatic, co- 
herent Langmuir waves were artificially injected as an initial conditio n instead of an electro n beam, 
yielding a strong signal of both fundamental and harmonic emission. lYu &: Guanglil (J2008) focused 
on the generation of backward-propagating waves through scattering (Eq[T]), confirming that di- 
rect creation of these waves is d ifficult in a particle-in-cell simulation. Finally, the simulations of 
J. I. Sakai fc M. Karlickvl (|2008l ) attempt to model the shock structure itself within the simulation 
box, in an attempt to represent the complete emission mechanism within one self-contained sim- 
ulation. Due to the inherent numerical limits given by a PiC code, both the accelerated particles 
and wave quantities are at odds with the other papers presented here. 
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Fig. 1. — Numerical simulation setup of the 2.5D PiC simulations. In addition to the homogeneous, 
thermal background plasma, two counterstreaming electron beam populations are present. 



Since the phenomenology of beam-driven wave instabilities depends on an electron distribution 
function that is far from equilibrium, MHD or fluid simulations are not suitable for this problem, 
but rather kinetic simulation methods have to be employed, with correspondingly larger computa- 
tional requirements. To obtain and analyze full spatial information about wave processes within 
the CME foreshock region (a s opposed to ID s atellite measurements), the kinetic particle-in-cell 
simulation code ACRONYM (jKilian et alj 120121 ) was used. This code, developed and maintained 
at the Department of Astronomy, University of Wiirzburg, is a fully relativistic, 2nd order particle- 
in-cell code for astronomical, heliospheric and laboratory plasmas. Using MPI-parallelization, the 
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code is running on all major supercomputing platforms. 



Due to the inherent length- and timescale requirements (governed by the plasma Debye length 
of about Xj} ~ 1cm) of kinetic simulations, it is impossible to model a complete CME, or even a 
significant part of the actual shock front within the bounds of a simulation. Rather, the microphysics 
of plasma wave interaction within the foreshock region is of central importance here. The actual 
simulations therefore model only the foreshock region, with the electron beam initialized as external 
input. 



Improving on earlier work (IGanse et al.ll2011bl ). the focus of this paper is the self-consistent 
creation of electromagnetic emission from Langmuir waves which have not been artificially created 
as monochromatic external input, but excited through instabilities within the simulation itself. 
Corre spondingly larger timescales and better spectral resolution are required over lKarlicky &: Barta 
(|2010|). The simulations are set up as 2.5D rectangular grids with periodic boundaries, which are 
homogeneously filled with the background foreshock plasma under quiescent coronal conditions 
(T ~ 1MK, p = 2.5 • 10 7 cm~ 3 , B = 1G), with a thermal particle distribution. On top of that, 
two counterstreaming beamed electron populations are added at t>Beam ~ 5f t h, wh ose density is 
10% — 15% of the total electron density (see Fig. [1]), based on in-situ ob servation data (Knock et al 



2001 



). Their pitch angle distribution is centered around 45°, following Karlicky Sz Vandas ( 20071 ). 



The selection of two counter streaming electron populations as oppos ed to a single b eam is 
motivated by our earlier results (jGanse et al.l l2011al ) , as well as work by lYu fe Guanglil (|2008l ) , 
showing that a creation of backward-propagating Langmuir waves of sufficient strength within a PiC 
simulation is itself numerically challenging. The counterstreaming beam setup provides symmetric 
creation of forward- and backward-pr opagating waves, and can physically occur in scenarios of 
multiple shock drift acceleration sites (jCairns et al.l l2003). 



To obtain physically valid results, the electron to proton mass ratio has been chosen as the 
physical m p /m e = 1836, the protons being part of the homogeneous, thermal background. In order 
to reduce the fundamental problem of numerical noise in our PiC-simulations, the particle number 
per cell had to be chosen quite large. Results have shown that below 100 particles per cell, nonlinear 
wave interactions were strongly disturbed by nonphysical effects due to statistical noise. This also 
sets a lower limit to the beam intensities which can be simulated given the current computational 
resources, as the generated wave signals have to be s ufficiently stronger t han the noise background. 
They are however still in reasonable agreement with lKnock et al.1 (|200ll ). 



Using large supercomputers, numerical simulation extents of 8192 x 4096 cells with 200 particles 
per cell were attainable, leading to physical extents of 80.28 x 40.14 metres. While this is still many 
orders of magnitude smaller than the expected size of the complete emission region, it resolves the 
kinetic plasma length scales and the wavelengths of all contributing waves with sufficient headroom. 

Multiple simulation runs were conducted to study the influence of multiple input parameters 
on the wave interaction behaviour in this scenario: The basic simulation setup outlined above was 
varied in beam density (a 50% increase) and background plasma density. 
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Fig. 2. — Energy distribution in two foreshock electron beam simulations. A transfer of energy 
into electric fields and subsequently into transverse magnetic fields is taking place, starting from 
t ~ 20c*;" 1 . In a direct comparison of two simulations with a 50% difference in electron beam 
density, a 52.7% increase in peak electric field intensity is observed. The resulting gain in magnetic 
field energy exceeds that, yielding a 172% increase. 



3. Results 



In the following, the results of two simulation runs, one with electron beam density 10% of the 
background plasma density, and one with 15% of background plasma density shall be presented and 
compared. Firstly, the spatially averaged distribution of energy into particle kinetic energy, electric 
and magnetic fields will be discussed, followed by decomposition of the field's energy contributions 
into individual spatial modes. 

The energy distribution within the simulations' constituents (Fig. [2]) yields that the initial 
energy, stored only in the particles' motion and the background B-field, starts to convert into 
electric field energy at around t = 20w~ . Soon afterwards, energy conversion into transverse 
B-field components is also observed. 

The initial transfer of energy into electric fields is attributed to the excitation of Langmuir 
waves by the presence of the electron beam, consistent with theoretical predictions. The peak energy 
content in th e electric field is ab o ut 2% of the initial electron beam's kinetic energy, confirming 
the results of iKarlickv &: Vandasl (|2007l ). The conversion of electric into magnetic energy (in the 
transverse components of the B-field) is then attributed to the creation of transverse electromagnetic 
waves. 
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A direct comparison between the simulation with the weaker electron beams (electron density 
10% nbg) and the stronger-beam simulation (ribeara ~ 15% n^g, same beam velocity), 
given in Table [3] shows that the increase in beam energy leads to a proportionate increase in 
peak electric field strength, and a much larger increase in transverse magnetic field energy. The 
increase in electric field energy by 53% is consistent with the linear relationship of beam strength 
to electrostatic mode excitation by Cherenkov-type instabilities. 



Beam electron density 
(as % of nb s ) 


Peak electric field energy 


Peak transverse 
mag. field energy 


10 


1082 erg 


0.18 erg 


15 


1652 erg 


0.49 erg 


Energy ratio 


1.53 


2.72 



Table 1: Comparison of two simulation runs, in terms of electron beam intensity and resulting 
energy content in electric and transverse magnetic field components. 

If the transverse magnetic field strength is actually caused by second-order wave couplings of 
electrostatic wavemodes, as outlined in section [H a quadratic increase in energy in this component 
relative to the electrostatic mode energy is expected: The observed increase by a factor of 2.72 
is in reasonable agreement with this prediction, as y/2.72 = 1.65 is comparable to the growth of 
the electrostatic mode by a factor of 1.53. Higher-order nonlinear effects might be the cause for 
additional energy excess here. 

In order to obtain more detailed quantitative information about the energy distribution into 
different wave modes within the foreshock plasma, the discretized simulation data was Fourier- 
transformed both in space, allowing analysis of instantaneous k-space spectra, and in time, to 
obtain dispersion diagrams. 

Figure [3] shows /c-space intensities of longitudinal E- and transverse -B-Fields at t = 65 uj~} . 
The x-axis corresponds to the electron streaming direction. 

By taking the Fourier transform in both temporal and spatial dimensions, the dispersion plots 
in Fig. [5] and [5] are obtained. They respectively show the relative intensities of the longitudinal 
electric and transverse magnetic fields in relation to u and k. 

3.1. Langmuir wave excitation 

The longitudinal E-field (Fig. [3l left picture) shows maxima of intensity along the streaming 
direction, as well as resonant excitations consistent with the streaming electrons' pitch angle dis- 
tribution, forming cones of high intensity around the x-axis. An additional, but weaker component 
perpendicular to the streaming direction is also visible. The cones around the x-axis are the ex- 
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Fig. 3. — /c-space distribution of wave energies in the foreshock plasma (arbitrary units). Left: 
Longitudinal electric field, Right: Transverse magnetic field. The x-axis corresponds to the electron 
beam direction. 



pected signature of beam-driven wave excitation, following the beam direction and being widened 
by the electrons' pitch-angle distribution. The components oblique and quasi-perpendicular to the 
electron beams' direction are not explained by thi s mechanism, and are as sumed to be scattering 
products of three- wave interactions, as outlined in lKarlicky &; Bartal (120101 ) . 

In Fig. [H k y = cuts through /c-space have been taken for both the weak beam and strong 
beam case, and Fourier-transformed in time to yield k — u dispersion diagrams, allowing for iden- 
tification of individual wave modes. In the right picture, intensity ratios between the weak beam 
and strong beam simulation are depicted, showing that the added wave intensity primarily leads to 
a broadening of the beam mode excitation towards lower frequencies. The dominant contribution 
of energy in these plots is located on the beam mode, with a constant phase velocity matching 
the average beam velocity (jWilles Cairnsll2000l ). Around the plasma frequency, where the beam- 
instability has its resonant maximum growth rate, the highest intensities are observed. The two 
symmetric, counter-streaming beams lead to symmetric dispersion behaviour. 

It should be noted that this beam mode is not an eigenmode of the background plasma, 
but rather presents an entropy wave transported by the beam component. As the beam weakens 
due to Landau damping, this mode has to undergo three-wave interactions to decay into actual 
plasma eigenmodes, as theoretically predicted in section [TJ Due to the short physical timescales of 
the kinetic simulations conducted here, a relaxation of the electron beams' bump on tail has not 
progressed to a large degree, so that the intensity in the beam mode is still dominant at the end 
of the runs. This stands in contrast to in-situ satellite observations, which are assumed to probe 
plasma conditions where most of the beam mode's energy has already transferred into Langmuir 
waves. Conversion efficiencies of beam to wave energies are therefore not directly comparable. 
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3.2. Electromagnetic wave emission 

The /c-space plot of the transverse magnetic field (Fig. O right) shows two separate emission 
regions (and their symmetrical counterparts): one at a low angle of about 20° against the beam 
axis, the other at near-perpendicular angles. Note that the signal-to-noise ratio of these signals, 
being produced by higher-order processes, is much lower than in the electric field components. 

From theoretical considerations of momentum and energy conservation ([T]), the near-perpendicular 
emission can be identified with the harmonic emission process, whereas the lower-angle emission is 
consistent with the fundamental frequency emission process. 



E x in x direction (longitudinal) Intensity ratios of longitudinal electric field modes. 
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k(1/cm) k(1/cm) 



Fig. 4. — Left: Dispersion plot of the longitudinal electric field intensity along the beam direction. 
The beam resonance deposits large amounts of energy into the beam driven mode. Right: Increase 
(ratio) of electric field spectral energy between the weak beam and strong beam simulation. 

The k x = plane of /c-space was used to obtain the k — lo dispersion plot in Fig. [5] (left). In 
the plot, the electromagnetic mode is clearly visible as a parabola with cutoff at the plasma fre- 
quency. Additionally, lower velocity modes dominate the low uj values, probably mainly consisting 
of improperly resolved Alfven waves (which are of no significance to the processes here). 

Most importantly however, two horizontal bands in this plot are noticeable, localized at u = oj pe 
and 2tOp e respectively. These bands do not follow the dispersion relations of any mode predicted by 
linear theory, which indicates that their origin is of non-linear nature. Kinematically, they correlate 
with the predicted three- wave interaction processes in Eq. [2] and [3] - assuming a coupling of waves 
of opposite-directed /c-vectors from the resonantly excited longitudinal modes. 

The /c-value of maximum coupling of the harmonic emission predicted by Eq. [8] is k y = 8.3 ■ 
10~ 3 cm, which matches the point of resonance visible in the simulation results (marked as a red 
arrow in Fig. [5]). It is noteworthy that the intensities of the two bands are quite similar, in contrast 
to observations which typically show higher intensity in the harmonic than fundamental emission. 
This is, however, assumed to be a propagation effect: due to the close proximity of the emission to 
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Fig. 5. — Left: Dispersion plot of the transverse magnetic field intensity (arb. units), with wave 
propagation direction perpendicular to the electron beam. Along the electromagnetic mode, reso- 
nances at uj pe and 2cj pe are visible (dashed lines). The theoretically predicted point of interaction 
is denoted by the red arrow. Right: Ratio of intensities between the weak beam and strong beam 
simulations. Significant increase in energy is visible at the fundamental and harmonic frequency 
bands, especially at their upper edge. 



the plasma frequency, the local medium has a much higher index of refraction for the fundamental 
than the harmonic emission. The harmonic emission can hence escape the emission region much 
easier, and leads to a stronger signal in observations. 

Apart from the plasma frequency and its harmonic, some intensity is also visible along the 
electromagnetic mode's dispersion relation at other frequencies. This excitation stems from the 
inherent numerical noise of the particle in cell simulation, caused by the relatively low particle 
number per cell. 



4. Conclusion 



Using our fully relativistic particle-in-cell code ACRONYM, we have simulated electron beam 
driven wave excitation and subsequent wave couplings in a CME-foreshock environment. The com- 
plete spatial and temporal information obtainable from a numerical simulation provides improved 
insight into the wave processes in the emission region over satellite measurements: Analysis of 
resulting k-space wave distributions yields kinematically sound emission patterns, consistent with 
analytical predictions about the type II emission mechanism. Further analysis of spatial and tem- 
poral data provides evidence of nonlinear three-wave interaction between beam driven electrostatic 
modes to create electromagnetic modes at fundamental and harmonic frequencies with roughly 
equivalent intensities. 
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The fundamental and harmonic em ission features vis i ble in /c-space representations of the 
simulated fields match those produced in lKarlickv Bartal (|2010l ). Since the emission in our runs 
was self-consistently produced from electron beams exciting Langmuir waves instead of artificial 
monochromatic excitation, these features are less strongly peaked, noisier and less pronounced, but 
otherwise share the emission characteristics. 



4.1. Outlook 

The beam intensity is just one of multiple tunable parameters of the emission process. Further 
parameters include beam velocity, background magnetic field intensity and the beam electrons' 
pitch angle distribution. Simulations are currently underway, studying the behaviour of the system 
under variation of these parameters. 

The enormous computational demands of particle-in-cell simulations are still a limiting factor, 
determining the maximum temporal extent of the runs. Slow wave interaction processes, like the 
scattering of Langmuir waves on sound waves (Eq. [1]) are still insufficiently resolvable, and future 
improvements in computing power and kinetic plasma simulation techniques will be required to 
create complete, all-encompassing simulations of this emission mechanism. 

An alternative simulation method worth studying would be a Vlasov code which, while equally 
or even more numerically challenging, would be able to simulate similar plasma behaviour free of 
effects due to statistical noise. 
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